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Abstract 

We consider simply generated trees, like rooted plane trees, and 
consider the problem of computing generating functions of so-called 
bare functionals, like the tree factorial, using B-series from Butcher's 
theory. We exhibit a special class of functionals from probability the- 
ory: the associated generating functions can be seen as limiting traces 
of product of semi-circular elements. 
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1 Introduction 

Let T n denote the set of rooted plane trees of size n. Simply generated 
trees are families of trees obtained by assigning weights uj(t) to the elements 
t € T = Un^n using a degree function ip(z) = 1 + Sfe>i 4 , kZ k (see [20]) . 
Basically, the weight u{t) of some t £ T is obtained by multiplying the 
factors rpfihi) over the nodes v of t, where d(v) denotes the outdegree of v. 
Our main topic is the study of generating functions 

associated with multiplicative functions B : T — ► R defined recursively by 
using a sequence of real numbers {-Bfc}fceN+ • We call such multiplicative 
functions bare Green functions: ^tgj" B(t)u(t) represents the sum of the 



Feynman amptitudes associated to the relevant diagrams of size n in some 
field theory, and the generating function is then a part of the perturbative 
expansion of the solution of some equation describing the system (see El 

EM )• 

In Section |1J we give an equation satisfied by Y when the weights 
come from some master function L(z) = ^2 m ^.Q L m z m , with = L(k)/k, 
\/k G N + . We use series indexed by trees, the so-called B-series, as defined 
in |13l I14j , to show in Theorem ^ that Y solves 

Y' = L(l + e)V(Y), 

where 6 is the differential operator 9 = zd/dz. 1 considers a similar prob- 
lem for additive tree Junctionals s(t) defined on varieties of increasing trees, 
like s(t) = 1xi(B(t)). Assuming some constraints on the degree function 
^f(z), it is proven that the exponential generating function 

S(z) = 5>(i) S (i)zKI/|i|!, 

is given by the formula 

S(z) = W'(z) [ (F'(u)/W'(u))du, 
Jo 

where F(u) = J2 m ^o M-^m) W m u m /m\ and W(z) = J2 m ^o W m z m /m\ solves 
W = ^>(yv). We also consider a central functional called the tree factorial, 
denoted by i! in the sequel, which is relevant in various fields, like algorith- 
mics ) stochastics jl 1 1 121) , numerical analysis (see for example ) , 
and physics We focus on its negative powers l/(t\) l+1 , I G N, which 

do not admit a master function when I ^ 1. p| solved the case I = 1 by 
using the so-called Butcher's group (see for example We provide 

in Theorem [2 a differential equation for the associated generating function, 
VI G N. 

In Section [3 we define special multiplicative functionals for which the 
weights Bk are related to the covariance function r of some gaussian pro- 
cess, as Bk = p 2 r(2k — 1), for some positive constant j3 > 0. We show 
that the generating function Y is related to the mean normalized trace of 
products of large symmetric random matrices having independent and in- 
dentically distributed versions of the process as entries. Theorem |31 gives 
then a differential equation for the evolution of the trace of a stationary 
Wigner processes. It follows that most of the examples given in [31 El can 
be expressed in terms of traces of large random matrices. In Section H3 we 
show how B-series can be useful for studying traces of triangular operators 
appearing in free probability. 
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2 Basic notions 



A rooted tree t G TZ is a triple £ = (r, V, .E) such that i) (V, E) is a non-empty 
directed tree with node set V and edge set E, ii) all edges are directed away 
from the root r £ V. The set of rooted trees of order n is denoted by TZ n , 
and the set of rooted trees is TZ = L) n 7Z n . A rooted plane tree t G T is a 
quadruple t = (r,V,E,L) satisfying i) and ii) and hi) L := {({w : vw G 
E}, L v ) : v G V} is a collection of \V\ linear orders. Given v G V, let 
ch(v) := {w : vw G E} be the set of children of v. d(v) := |ch(i>)| is the 
outdegree of v. A rooted planar tree can be seen in the plane with the root 
in the lowest position, such that the orders L v coincide with the left-right 
order. Next consider the partial ordering (V, ^) defined by u ^ v if and 
only if u lies on the path linking r and v. Given v G V and t G TZ let 
t v be the subtree of t rooted at v spanned by the subset {w; v ^ w}. A 
rooted labelled tree is a quadruple t = (r,V,E,l) satisfying i) and ii), with 
a labelling I: V\{r} — ► [\V\] := {1, • • • , \V\} such that l(u) < l(v) when 
u < v. The set of rooted labelled trees of order n is denoted by C n . Let 
C = U n C n . This family is a special variety of increasing trees, as defined in 

mm. 

We next assign weights to the elements of T n , the set of rooted planar 
trees of order n: the resulting family of trees is said to be simply generated 
(see (25)- Given a sequence i\) = {tpk}k€N of real numbers with tpQ = 1, 
define recursively the weight uj{t) of t G T as 

k 

i=l uGV 

where t±,--- ,tf- are the d(r) subtrees of t rooted at ch(r). Let ip(z) := 
1 + J^'kLi ^k^ be the generating function of the weight sequence ip. Our 
favourite example is ip{z) = 1/(1 — z), with u(t) = 1, Vi (see ^ 120] for 
various interesting choices). 

We will be concerned with functionals B : T — > M., where T = U n J- n , 
called bare Green functions. This terminology is taken from quantum field 
theory where bare Green functions occur during the action of the renormal- 
ization group (see for example [HJ, § 4.2 or , § 6.1). Let B denote the set 
of bare Green functions. Any element B G B is given through a sequence 
of functions B^ : R — > R, k G N + , which are usually Laurent series in 
some variable x (see for example 8J. In what follows, we simply write the 
sequence as {-Bfc}fceN+- 

Definition 1 The bare Green function B G B, B : T — > K, associated 
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with the sequence of functions {-Bfc}fceN + is defined recursively as 

k 

£(t)=Bi t in B (*i)> 

i=i 

where ti, ■ ■ ■ , ifc are the d{r) subtrees oft rooted at ch(r), and where \t\ 
denotes the number of nodes oft. 

Notice that the value of B at t £ T does not depend on the linear orders 
and is independent of the labellings. When dealing with rooted trees, we 
will adopt the notation t = B + (t\, ■ ■ ■ , t^) for the operation of grafting the 
rooted trees t\,--- ,tk, that is by considering the tree t obtained by the 
creation of a new node r (the root) and then joining the roots of t\, ■ • ■ ,tk 
to r. Bare Green functions appeared also in the probabilistic literature in 
specific situations. The basic example, in algorithmics in numerical 

analysis (see [U E] ) , in stochastics EOj and in physics (see for example 
[0] ) is the tree factorial, defined by 

Definition 2 Let t £ TZ with t = B+(t\,--- , Then the tree factorial 
t! is the functional B £ B defined by t\ = |£| IliLi associated with the 
sequence {Bk} given by B^ = k. 

Remark 3 It should be pointed out that the functional acting on trees, given 
as s(t) = \n(B(t)), for B £ B with B^ > 0, Vfc ^ 1, is an inductive map or 
an additive tree functional, as defined in Interestingly, B(t) = 1/tl is 
used in UMj to define a probability measure on random search binary trees, 
and [77^ provide precise asymptotics for ln(i!). 

3 Generating functions 

We first give some basic results on tree factorials, symmetry factors, and 
generating functions associated with bare Green functions. 

Definition 4 Let t £ TZ. Then a(t) is the number of rooted labelled trees 
t' £ C of shape t £ TZ, where the shape of a labelled tree (r, V, E, I) is (r, V, E), 
n(t) is the number of rooted plane trees of shape t, and o~(t) is the symmetry 
factor of the tree, to be defined later. Moreover, let loc be the weight function 
associated with elements of C, with weights given by ipp. = l/k\. 

Notice that a(t) is the Connes-Moscovici weight in quantum field theory 
(see [Sill!)- The symmetry factor satisfies the recursive definition: 

a({r}) = 1, 
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a(B + (t^, ■ ■■,%»)) = nMh) ni ■ ■ ■ n k \<j{t k ) n \ 

where the indices rtj means that t is obtained by grafting n\ times the tree 
t\, and so on, where we assume that the ti are all different as rooted trees. 

Lemma 5 Let t £1Z. Then 

<*{t)*(t) = ^, (1) 

and 

a(t)t\ = \t\\uc{t)n{t). (2) 

Proof: (Pi is well known (see for example 0). Suppose that t G 7£ is such 
that i = B + (t T l 1 , ■ ■ ■ ,i^ k ), the grafting of ni times the tree ti, and so on, 
where we set that the trees t±, ■ ■ ■ ,t k are different as rooted trees. Then 



<t) = (ni V + ? )! ^r ■ ■ ■ «t k r> 

ml •••rafc' 



Using the recursive definition of u)(t) and the definition of u>c, we have 

(m H h ra fe )! 



Therefore 



uc{t)n{t) ' ojc(h)K(ti) U£(tk)K(tk)' 

and the results follows from the recursive definition of the symmetry factor. 

□ 

Then 

£ B(t)u(t) = £ (3) 
' — ' ' — ' OJr \t\\ 

where we have used (J2J of Lemma |5J 
Consider the generating function 

Y ^ = E ^ E a(*)S(t)tIw(t)/ WjC (t). (4) 

Given t £ 1Z, the ratio lo/ujc is associated with the weight sequence = 
V'fcA;!; using the expansion ifj(z) = 1 + Y^k^i = 1 + Ylk^i^k/k\)z k , we 
see that ^ = ?/^ fc )(0). Consider the elementary differentials 5 (see Section 
EJ defined by 
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Definition 6 

k 

5 W = 1, 8 t = ^ k \0)T\5 u , —=S, 
T=i ^ 

when t = B + (ti,--- , ij), where * denotes the tree of a single node. For 
a map a : TZ U {0} — ► R, a formal power series of the form Y(z) = 
a(0)yo + Ete^" z '*' a ( t )^ a W/N ! * s ca/Zed a B-series H31 [7^] /. 

Remark 7 When B(t) = t\, the series Y is given by 

Y{z) = Y J {"{t)luc{t))z\ t \j\t\\. 
tec 

Setcpk = ipkkl, Vfc, and consider the degree function <fi(z) = ^+Yl,k>i^ ( l ) k/k\)z k 
Following [1], Y solves Y' = (f>(Y) (see also [20]). We shall see in the next 
section that it is a natural consequence of B-series expansions of solutions 
of ordinary differential equations. 



4 Runge-Kutta methods for functionals over trees 

Consider a dynamical system on R 

^X(s) = F(X(s)), X(s ) = X , 

for some smooth F : R — * R. The solution of this dynamical system has a 
B-series expansion of the form 

X(s) =X + ^2 ( 8 ~*f® a(t)6t(8 ) t 
tell ' 

where the elementary differentials 5 is defined recursively by 

d k F 

o~{r} = f{so), d~t = -gkjK ■■■S tk , 

when t = B + (t\, ■ ■ ■ ,tf.). These kinds of expansions have been treated in 
great detail in [I] and [2] and developped independently in combinatorics 
(see for example El ) • Suppose that sq = for simplicity. Butcher con- 
sidered what happens with numerical approximations of the exact solution, 
the Runge-Kutta methods, which are themselves B-series |13l I14j : here we 
focus on specific B-series, which are associated to bare Green functions. Let 
B £ B be such that there exists a power series 

L(Z) = L mZ m , 
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with 

S t = ^, Vfc6N+. (5) 
k 

Bare Green functions satisfying (JSJ are used in practical situations in quan- 
tum field theory (see 0, § 4 and |Hj , § 6.1). Consider Euler's operator 
6 = z(d/dz), with P(9)(z n ) = P(n)z n , Vn G N, for any polynomial P, and 
consider the formal operator L{6 + 1) acting on monomials as 

L(9+l)(z n ) = J2 L ™(9 + l) m (z n )= ^i m (n + l) m z™ 

m^O m^O 

= L(n + l)z n . 

Given a power series Y(z) = ^2 irri> ^a m z m converging for \z\ ^ 1, we can 
define L(8 + l)(Y)(z) := X^m>o a m,L(m + l)z m , which converges for \z\ ^ 1 
when the sequence (L(k))k^i grows sub exponentially. We will not focus 
on convergence questions here, and work at the formal level. Let B be a 
bare Green function with weights (-Bfc)fc>i> such that (JSJ) holds for some 
power series L. It should be pointed out that deal with the master 

function L, but do not give explicitely an equation for Y. The next Theorem 
provides an equation; its proof uses explicitely B-series. 

Theorem 1 The formal power series 

n^^i^B^, (6) 

tell ' 

solves Y' = L(l + 9)ip(Y). 
Proof: 



fc^O ' (t v ..t k )£-Rk 1 11 • 1 i=l 

For given (ii • • • t k ) £ TZ k , set t = B + (h,--- ,t k ). Then 1**1 = 1*1 ~ ^ 
V> {fc) ( )<^i ■~ s t k = 5 ti and S (*0 " • •- B (*fc) = -B(*)/J3|*i- The associated term 
becomes 

|t|-i-B(*h a(ti)-afe) *' 
Z * 1*1' 

Next, every rooted tree t € TZ can be decomposed uniquely as t = B + 
(i™ 1 , • • • , im m )> meaning that i is obtained by grafting n\ times t\ and so on, 
where the U are different as rooted trees, with k = n% + ■ ■ ■ +n m . Collecting 
the terms associated with t, we get the contribution 

z^Bjt) tl A a(t' 1 )...g(t' k ) 

k\ Bu\ *|t| ^ l*i|!---|*i|! ' 

1 1 (t' 1 -t' k )ell k 111 1 kl 
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where * means that the sum is taken over all the collections (t^ ■ ■ ■ t' k ) £ 7Z k 
such that t = B + (t' 1} ■ ■ ■ ,t' k ). The above sum reduced then to 

(m + • • • + n m )\ a(h)^ ■ ■ ■ a(t m ) n ™ _ 1 ( 1 )Wl _ 1 NTlm 



ni\---n m \ k\\ti\\ ni ■ ■ ■ \t m \l nm m! a{tx)t\\ n m \ a(t m )t m \ 

where we use the first identity of Lemma OD Using the recursive definition 
of the symmetry factor a, we obtain 



E 



a(<i) ■■■a(t' k ) 



1*1 1 !*!«(*) 



t\ a(t) \t\\ ' 

We thus get that the contribution associated with t £ 1Z is given by 
*! B|t| 1*1 l<|! B m • 1 

Therefore 



1*1-1 



teft I 1 ! 1 1 

£(^rT)T B < f > (!5 <=d7' 



□ 



Remark 8 ^4s we Ziawe observed in Remark\^ s(t) = 1n(B(t)) is an induc- 
tive map when the weights are positive. It turns out that the exponential 
generating function associated with s can be given as an integral transform 
for varieties of increasing trees (see for example Section^). This is the topic 
ofW- 

Example 9 

WhenL(z) = z, with B^ = 1, and^(^) = 1/(1 — z), one has YlteTn+i = 
C n , the Catalan number of order n, with C n = ( 2n ) /{n+ 1). Then Y(z) = 
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z X^n>o z "Cn is solution of the differential equation Y'(z) = L(l + 0)(1/(1 — 
Y{z)f), that is Y'(z) = (z/(l - Y(z)))'. The unique solution with Y(0) = 
satisfies Y(z) = z/(l — Y(z)), or Y(z) = (1 — y/1 — 4z)/2, corresponding to 
a well known result. 

□ 

[5] , § 5.3, considers the case where B(t) = (1/tl) 2 , which is not of the form 
given in ©: in this situation, B k = 1/A; 2 , with L(z) = l/z. The solution is 
obtained by using the stucture of the so-called Butcher's group of B-series 
(that is series of the form © , where the group structure in given in ^1 ) 
by tensoring known B-series: 



Example 10 

Consider the bare functional given by B k = 1/k 2 , with B(t) = l/tl 2 . Fol- 
lowing Brouder, the associated B-series, as given in @, is solution of the 
second order differential equation 

Z Y" + Y' = ip(Y). 

When ifj(z) = exp(z), the solution is given by 

00 n i 

Y(z) = -2ln(l-z/2) = J2 Z -^z- 1 , 

n=l 



giving 



^ a{t) _ (n- 1)! 



□ 



We study the general moment problem B(t) = (l/tl) l+1 , I G N, by 
working directly on a suitable differential equation as follows: the operator 
L{9 + 1) takes the form L(6 + 1) = 1/(9 + 1)'. Assume that the differential 
operator L(9 + 1) is invertible. Then the formal systems becomes 

L(9 + l)- 1 ^Y = i>{Y). (7) 

Consider again the second moment problem for the tree factorial, with B k = 
1/k 2 and L(k) = 1/k. Choose L such that L(z) = l/z; the inverse operator 
might be equal to L(9 + =0+1, and, if this is the case, 

(e + l)Ay = ^(y), 
with (0 + l)(d/dz) = z(d 2 /dz 2 ) + (d/dz), see Example ITU1 
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More generally, if one considers the moment of order I + 1 E N of the 
inverse tree factorial, the choice L(k) = l/k l should give (6 + l) l -^Y(z) = 
ip{Y)- Our result, Theorem [2] below shows that the formalism of inversion 
is correct in term of power series. This result sheds some light and extends 
the computations done in J2] for the second moment, and its proof avoids 
computations in the Butcher's group. 



Theorem 2 The B-series Y(z) associated with the moment of order (I + 1) 
of the inverse tree factorial satisfies the differential equation 

(e + l) l ±Y = ^(Y). (8) 

Proof: Let 

t&TZ 1 1 

Then 

™ = ^-T 2. | tl |!...| tfc |! ( tl !...t fc !)' 

For given {tf-t k ) G K k , set i = ,t k ). Then £\ W = 1*1 - 

1, ip^(0)5 tl ---5t k = S t , and (tx\---t k \) 1 = t\ l /\t\ l . The associated term 
becomes 

z^\t\[ ajh) ■ ■ ■ a(t k ) 
tV 1 |ti|!...|t fc |! ' 
Proceeding as in the proof of Theorem Q we get that the contribution asso- 
ciated with i G 7£ is given by 

z^- x \t\ l . \t\a(t) z^- 1 \t\ l . 

— in— *^jr = w^.w am - 

On the other hand, 

{e+l)l iz Y{z) = (^^TrtW^ 



^iJ^ltl-i i 



□ 

In the next section, we show that traces of certain products of Wigner 
matrices (see for example 24,) provide natural examples of bare Green func- 
tions. 
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5 Wigner processes 

Definition 11 The N -dimensional random matrices Tjv := {'Ji,j)i^i l j^.N 
are called Wigner matrices of variance j3 2 if the following holds. 

• Each Tn is symmetric, that is, jij = 7^. 

• For % ^ j, the random variables jij are independent and centered. 
. Fori +3, E( 1 2 j )=p 2 . 

• For any k ^ 2, E(|7jj| fc ) ^ c^, where cj~ is independent of i ^ j. 

Definition 12 The sequence T^(k) := (7ij(fc))i^jj^7v °f N -dimensional 
random matrices, indexed by k ^ 1, is called a Wigner process of variance ft 2 
and correlation function r, r(k,k) = 1, \r(k,m)\ ^ 1 and r(k,m) = r(m,k) 
if the following holds. 

• Each Tjv(fc) is a Wigner matrix of variance ft 2 in the sense of defini- 
tion \tl\ 

• For i ^ j, each process (7i,j(&))fc is independent of the others. 

• For i 7^ j, the process ( 7?, j (&))/<; is r -correlated, that is, for any k ^ m, 

E(7ij(A;)7ijM) := f3 2 r(k,m). (9) 
A Wigner process is stationary when r is such that r(k,m) = r(\k — m\). 

Let -Djv be a sequence of random diagonal matrices, with independent and 
identically distributed entries of law fi, having finite moment [i^ = /j(X k ), 
k ^ 1, with /ii = 1. Let 

k 

: =N- k / 2 D N H(T N (m)D N ), 

771=1 

and set 

5f(r)=iV- 1 E(tr(Qf)). 



Involutions, Dyck paths and rooted plane trees 

For k ^ 1, [k] := {1, 2, . . . , k}, I(k) is the set of the involutions of [k] 
with no fixed point, J{k) is the subset of T{k) of the involutions a with 
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no crossing. This means that the configurations i < j < a(i) < a(j) do 
not appear in a G J(k). Let % G cr(cr) denote the fact that i < a(i). Let 
D(2k) be the set of the Dyck paths of length 2k, that is, of the sequences 
c := (cn)o<n<2fc of nonnegative integers such that c = c 2 k = 0, c n - c n _i = 
±1, n G [2k]. Thus, exactly k indices n G [2k] correspond to ascending 
steps (c n _i,c n ), that is, to steps when c n = c n _i + 1. We denote this by 
n G asc(c). The k others indices correspond to descending steps, that is, to 
steps when c n = c n _i — 1, and we denote this by n G desc(c). We make use 
of bijections between V(2k) and J (2k) 2 . If c G V(2k), <j)(c) := a G J{2k) 
is an involution which maps each element of desc(c) to a smaller element 
of asc(c). Thus, cr(<r) = asc(c). More specifically, if n G desc(c), cr(n) 
is the greatest vn ^ n such that (c m _i,c m ) = (c n ,c n — i ). Finally, the set 
V(2k) is in bijection with J^k+i-, the set of rooted plane trees on k + 1 nodes, 
where the bijection is given by the walk on the tree from the right to the 
left (see for example ;25J. Let at denote the involution of ^7(2(^1 — 1)) 
corresponding to t G T . Given t G J~k+i, consider the walk on t from the 
right to the left: every edge (v, w) with w G ch(t>), is crossed at some instant 
s v G [2k] as (v — > w) and at a later time s w G [2A;] as (w — ► v). Clearly, 
s w = s v + 2(1^1 — 1) + 1, where t w is the subtree of t rooted at node w, that 
is the subgraph of t induced by the nodes u with u ^ w. at is such that 
&t(sv) = Sw and vice versa. 



Fig 1. Bijections between J-k+x, T>(2k) and J(2k) (s v = 2 and s w = 5). 

Proposition 13 Assume that the covariance r is such that r(l,m) = r(\l — 
m\). Then, the functional B r G B given by the weights 




w 




2 4 6 8 10 12 




B r k = (3 2 r{2k-l), Vk ^ 1, 



is such that 



B r (t) 



J] {(5 2 r{i,a t (i)))- 




iEcr(at) 
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Proof: Let t G T . Let s v < s w be the instants where the oriented edges 
(f — > w) and (w —> v), w £ ch(v), are crossed during the walk on the tree. 
r(s v ,a t (s v )) = r(s w - s v ) = r(2\t w \ - 1), and thus (3 2 r(s v , cr t (s v )) = B^. 

Finally, Uiecr(* t ) P 2 r{i,Q- t {i)) = Enroot B \t w \ = Br (t)/ B \t\> as required. 

□ 

As we have just seen, every Wigner process with covariance r such that 
r(l,m) = r(\l — m\) produces a bare Green function B r G B. The converse is 
not true, that is, there exists B^B such that B is not of the form B = B r for 
some covariance function r. Set B w = {B G B; 3 a covariance r with B = 
B r }. 

Let Vv be the generating function of the weight sequence i\) k = Mfc+ij 
and let t G T be the associated weight function. 

Theorem 3 Let (Tiv(fc))fe>i &e a stationary Wigner process of covariance 
function r and variance j3 2 , and let be a sequence of random diagonal 
matrices, independent of the Wigner process, with i.i.d. entries Xj of law \x, 
with hi = //(A) = 1 and finite moments = )> Vfc. Then 

B? k (r) — B 2k (r) = -L_ ]T ^(t)^(t), 

and -B2fe+i( r ) — * 0, -/V —> oo. Assume that the covariance is such that there 
exists a power series L r {z) with B T k = L r (k)/k, \/k. Then the formal power 
series 

Y{z) = Y J ^BlB 2{k _i ) (r), 

solves 

r = u\e + \)^{Y). 

Moreover 

J2z k B 2(k _i ) (r) = z^(Y). (10) 

Example 14 

Let B(t) = l/t\. If a tree t has n nodes and n — 1 edges, then the re- 
quirement B n = 1/n is satisfied iff /3 2 r(2n — 1) = 1/n, that is r must be 
such that P 2 r(k) = 2/{k + 1), k G 2N + 1. By construction, r(0) = 1 
and therefore (3 2 = 2. l/(x + 1) is positive definite, which implies that 
B(t) = l/t\ is element of B 1 ". Next, from Theorem^ the generating function 
Y l z ) = Ete^'^W*) is solution of the system (d/dz)F(z) = 
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Assume that \i is the point mass Si, that is each matrix Dn is the iden- 
tity matrix of size N, with t/Ju,(z) = 1/(1 — z). The solution of the sys- 
tem is Y(z) = 1 - VI - 2z = 2Y(z/2), where y is the series given m 
Example |HJ On the other hand, Proposition and Theorem |3] show that 
Y{z) = X)jfc>i zk ^k^2(k-i)( r )- Therefore the limiting mean normalized trace 
B2k{r) of the product of correlated random matrices Y\m=i^N{fn) is such 
that B2k{r) = E(Z 2fc )//c!, where Z denotes a normal N(0,1) random variable. 

□ 

Example 15 

Consider as in Example |H1 the special case where L{z) = z. The associated 
inductive parameter (see Remark [SJ is the tree size. The covariance r is 
constant with r(k) = 1, and Bl = 1. Then the generating function Y is 
solution of the fixed point equation Y(z) = zip^(Y(z)) (either by Theorem 
H or by (|l()jl). Notice that in this situation, Tn(iji) = rjv(l), and thus 
Bu (r) describes the mean normalized moment of the spectral measure of 
the random matrix Dj^(Tj^(l)Dj^) k . This example can be extended by 
considering L(z) = zp z , for some < p ^ 1. When Djy is the identity 
matrix, Y(z) is related to the Rogers- Ramanuj an continued fraction |19j . 
and corresponds to the generating function associated with path length, see 

pug. 

□ 

Proof of Theorem [3J- The first part is a generalization of Theorem 1 of 
[T9] , Set 7ij(m) = jij(m)Xj, and IV (m) = T^(m)Dpf. The mean nor- 
malized trace adds the contributions E(i) = E(Aj 7j j 1 • • • ; ji k _ 1 i k ), for paths 
i = (ii)o^i^ki with i\ G [N] and io = i^. The jij are centered, so that any 
edge appearing once appears at least twice. Given i, define £i = 1 and 
ei = +1 when ii g" {^q,-- - , H-i}, and e; = — 1 otherwise, and consider the 
walk c = (q) defined by q = Y^j=i £ j-> with Ck ^ 0. The support of i is 
s(i) = {if, ^ Z ^ A;}, of size s = \s(i)\, with s ^ 1 + The contribu- 
tion E(i) is independent of the labels if, they are N(N — 1) • • • (JV — s + 1) 
labellings giving the same walk c, with the same contribution. Thus, the 
normalization _/V~( 1+fe / 2 ' shows that the only walks surviving in the large 
N limit are those with s = 1 + k/2. This shows that B^ (r) — ► when 
k is odd. Concerning B^r), s = 1 + k means that every edge occur- 
ing in the path occurs exactly twice, in opposite directions, c is a Dyck 
path of V(2k); let t G T be the associated rooted plane tree, with in- 
volution at. Using the right to left walk on t and the independence of 
the random variables, the contribution E(i) of any path leading to c or t 

is E(i) = n mec , (CTt ) E (7(^)7(^(m)))E(n,A^ )+1 ) where d(v) = \ch(y)\. 
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From Proposition El one obtains E(i) = (B r (t)/B^. + x ) T[ v with 
B 2k(r) = Y,ter k+ S Br ^/ B l+i)Y\-v^d(v)+i, as required. (JTHJ) is a conse- 
quence of the multiplicative form of bare Green functions and of Lemma 
1.9, chap. III.l of H3|. 

□ 

These results show that the elements of appear naturally in the 
computation of normalized traces of products of large random matrices (see 
for example [23 ). In the next Section we illustrate B-series by considering 
triangular operators from free probability. 

6 On Dykema-Haagerup triangular operator 

Let B be an algebra and A be a B bi-module. Let k : A x A — > B be a 
bilinear map. We follow [22] by defining the product a\ » K a 2 = K(a\,a 2 ), 
ai,a 2 E A, and setting 

i) (bai) m K a 2 = b{a x m K a 2 ), 

ii) (aib) » K a 2 = a\ » K (ba 2 ), 

Hi) at m K (a 2 b) = (ai m K a 2 )b. 

Let a E J{2n) be an involution of [2n] without fixed point and without 
crossing. Given a word a = a\ ■ ■ ■ a 2n in A, a induces parentheses on a, 
and the preceedings rules permit the evaluation of this parenthized word. 
This extends to a map K a on A 2n . Sniady defines such maps to prove a 
conjecture of Dykema and Haagerup on generalized circular elements. Let 
(B C A, E) be an operator valued probability space, that is A is a unital 
*-algebra, Bcian unital *-subalgebra and E : A — * B be a conditional 
expectation (linear, E(l) = 1, and E(biab 2 ) = biE(a)b 2 , V&i, b 2 £ B, a E A). 

Definition 16 T E A is a generalized circular element if there is a bilinear 
map k satisfying the rules i), ii) and Hi) such that 

E(b 1 T s 'b 2 T S2 ---b 2n T^)= K*Q>xT",.-- ,b 2n T s ^), 

aej(2n) 

E^T^&sT 82 • • • b 2 n+iT S2n+1 ) = 0, 
V&i, • • • , b 2n+1 E B and Vsi, • • • , s 2n +i E {1, *}. 

The triangular operator T of Dykema and Haagerup is obtained from B = 
C[x], the *-algebra of complex polynomials of one variable by setting 

[ K (T,bT*)](x) = f b(s)ds, 
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[k(T*,VT)](x) = / b(s)ds, 
Jo 

[K(T,bT)](x) = [K(T*,bT*)](x)=Q. 



T is the limit for the convegence of *-moments of large upper triangular 
random matrices Tjv ( jlUQ- Define a trace r as (see \22\ ) 

r(a) = r(E(o)), r(6) = / 6(«)da. 

JO 

In what follows, we use P-series (where P stands for partitioned dif- 
ferential systems, see |H3). We follow jS], and adapt his notations to 
P-series. Given some function ip, and two kernels {a x (u,v)) uv& \Q t x] and 
(a y (u, u)) U) „ € r 0j i] , consider the iterated integrals ^ and $i which are func- 
tional over 1Z defined by = </>u(*) = L and, for i = B + (t±, • • • , 



Z(t)=f[ [ arfavWMdv, 



k J 



^(*)=II / a"(«,«)^(*i)dt;. 

Lemma 17 Xe£ a x (u,v) = I[o, u ](i>) arad a y (u,v) = Ir w u(«). Then 
r(TTT= E jf = E jf 



Proof: The word = (TT*) ■ ■ ■ (TT*) is of the generic form with 61 = 
• • • ^2n = 1 ( Definition I16|) . Let i G ^"n+i with associated involution o~t (see 
Section |SJ). Let s v and s m be the instants where the walk on t crosses the 
oriented edges (v — > w) and (w — ► v), with w G ch(u). We colour these edges 
by giving colour '1' to (v — ► w) when the symbol in W located at position 
s v is T, and give the colour '*' otherwise. Clearly, both edges have different 
colours, and the elements of the set of edges {(v — > w); w G ch(u)} (the 
children of v in t) have the same colour. The result is then a consequence 
of the definition of the product with the rules i), ii) and iii). 

□ 



Remark 18 Iterated integrals are naturel objects to consider in the setting 
of Butcher's Theory. For example, in the framework of Theorem the 
iterated integrals <j) u {t) defined by 4> u {t) = Y\ i=i Jq L(0 + l)((p v (ti))dv, when 
t = B + (t±, • • • , tk), are such that (f>i(t) = B{t), Vi G T . 
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Proposition 19 The P-series 

X u (s) = X + Y,—a(t)t\5 t / ^(u^^di;, 
ten ' '' ^° 

and ^ ^ 

ten ' ^° 
are solutions of the integral system 

X u {s)=X + s [ a x (u,v)t/j(Y v (s))dv, 
Jo 

Y u {s) = Yx + s f a y (u,v)^(X v {s))dv. 
Jo 

Proof: This is consequence of Butcher's general theory (see (!]). To prove it 
more directly, proceed as in the proof of Theorem ^ 

□ 

Corollary 20 Let Xq = Y± = 0. Assume that a x (u,v) = I[o u i(w) and 
a y (u,v) = Ir U)1 i(u). Suppose that ifi(z) = 1/(1 — z). Then 




-i 



This result shows that the generating function of the *-moments of the 
operator TT* can be obtained by solving the system given in Proposition 
1191 We recover in this way a result of |1U) . Lemmas 8.5 and 8.8. 

Lemma 21 In the setting of Corollaru \2(A the generating function Yq(s) 
solves 

G( - S VT ^ ) = s, (11) 
1 - Yo(s) 

where G(z) = zexp(— z), that is, L(s) = s/(l — Yq(s)) and G are inverse 
with respect to composition. Moreover T(TT*) n = n n /(n + 1)!. 

Proof: We solve the integral system by looking for solutions of the form 
X u (s) = 1 — exp(Au) and Y u (s) = 1 — exp(A'(u — 1)), with (d/du)X u (s) = 
s/(l - Y u (s)) and (d/du)Y u (s) = -s/(l - X u (s)). We deduce that A' = -A 
is solution of the equation A + s exp(— A) = 0. The formula for the moments 
of TT* is a consequence of Lagrange's inversion formula. 

□ 
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